Nonlinear von Neumann-type equations: Darboux invariance and spectra 



Maciej Kuna 1 , Marek Czachor 1 ' 2 , and Sergiej B. Lcblc 1 
1 Wydzial Fizyki Technicznej i Matematyki Stosowanej 
Pohtechnika Gdanska, id. Narutowicza 11/12, 80-952 Gdansk, Poland 

2 Arnold Sornmerfeld Institute for Mathematical Physics 
Technical University of Clausthal, 38678 Clausthal-Zellerfeld, Germany 



0\ 



S3 



> : 

m . 

(N ■ 

o ■ 
o : 

00 ■ 

^ ■ 
Ok 

> ■ 

cr. 



Generalized Euler-Arnold-von Neumann density matrix 
equations can be solved by a binary Darboux transformation, 
given here in a new form: p[l] = e PhlW ^ pe~ Plll ^ /v \ where 
P — P 2 is explicitly constructed in terms of conjugated Lax 
pairs, and p, v are complex. As a result spectra of p and p[l] 
are identical. Transformations allowing to shift and rescale 
spectrum of a solution are introduced, and a class of station- 
ary seed solutions is discussed. 



I. INTRODUCTION 

In this Letter we want to address the following prob- 
lem. Assume we have a nonlinear von Neumann equation 
(vNE) whose solutions can be generated by a binary Dar- 
boux transformation [jl]-[4j. The question is what is the 
relation between spectrum of a "seed" solution p and this 
of the Darboux-transformed p[l]. The issue is especially 
important if p[l] is to be a density matrix: Its spectrum 
cannot be negative, but does the fact that we start with 
a positive p imply positivity of p[l]? 

Quite recently density matrix solutions of the nonlin- 
ear Euler-Arnold-von Neumann equation (EAvNE) ip — 
[H, p 2 ] were constructed by this technique in || . The 
solutions p[l] were obtained from the seed p satisfying 
[A a ,H] = 0, a £ R, with nontrivial A a := p 2 — ap. The 
explicit form was shown to be 



/O [l](t) = e - iaH *(p(0) + ( M -/i)^(f)- 1 
[\<p(0))(<P<P)\,H 



ei A a t\ e iaHt 



(1) 



where p(0) is an initial condition for p, p is a complex pa- 
rameter of the Darboux transformation, \(f(0)} a solution 
of the Lax pair at t = 0, and 

F a (t) = (<p(Q)\exp (ifj^Aj) \<p(0)}. 

An analysis of explicit examples (e.g., for H correspond- 
ing to a 1-dimensional harmonic oscillator) revealed sev- 
eral interesting and unexpected properties of p[l](t). For 
example, in all the examples the dynamics was found to 
be asymptotically linear, but with different asymptotics 
for t — > +oo and t — > — oo and a kind of "phase transi- 
tion" around t — 0. Also in all the examples spectrum 
of was equal to this of p(0). It is difficult to see 

these properties in the formula (|l|) and therefore it was 
not clear to what extent they followed from the choice of 



the seed p, and whether they are a feature of only this 
particular nonlinear equation. 

In this Letter we prove that spectra of p[l] and p coin- 
cide whenever the equation they satisfy can be obtained 
as consistency condition for a Lax pair which is covariant 
under the binary Darboux transformation. The result is 
therefore not a pecularity of the EAvNE and does not 
depend on the choice of seed solutions. We shall give a 
general proof which does not refer to any concrete equa- 
tion but the notion of the binary-Darboux covariance will 
be illustrated on a class on nonlinear vNE's. Finally we 
shall show how to shift and rescale the spectra. Such 
properties are important in order to generate normalized 
density matrices p[l] from p's that are either nonpositivc 
or nonnormalizcd. 



II. ZAKHAROV-SHABAT PROBLEM AND 
DARBOUX- TYPE TRANSFORMATIONS 

Let V and J be linear operators, s a parameter, and 
p a complex number. We shall consider two types of 
Zakharov-Shabat (ZS) problems. Let \<p) G TL be a vector 
from a Hilbert space. The Hilbert-space ZS equation is 



id s \y) = (V-pLJ)\<p). 



(2) 



The second type of equation is obtained by taking a 
"second-quantized" problem 



id s if = (V — pJ)(p 



(3) 



where (p is an operator. 

In order to introduce the binary Darboux transforma- 
tion of a ZS problem we supplement (||) or (||) by an 
appropriate pair of conjugate ZS equations: 



id s (x\ 
-id s (i/j\ 



<x\{V-vJ), 

W(v-xj), 



for (|), and 



- id s x = x(V - vJ), 
-id s yb = rjj(V - A J) 

for (||). We define the operator P by 

p _ hgXxj 
(x\<p) 



(4) 
(5) 

(6) 
(7) 



(8) 
(9) 
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for the two cases, respectively, where the operator p 
satisfies p 2 — p and d s p — 0. The inverse in (||) is 
understood as the one in a p-invariant subspace: For 
any operator x — pxp its inverse x^ 1 — px~ 1 p satis- 
fies xx^ 1 = = p. In both cases P is an idempotent. 
For example, 

p 2 = 0(px<pp)~ 1 x ( p(px ( pp)~ 1 x 

= ( p(px ( ppy 1 px ( pp(px'fpy 1 x = tpp(px ( pp)~ 1 x = p- 

Define 

V[1]=V + {p-v)[P,J] 



and 



V — fl 



A 

v - 



P 



A — fi 



P 



(10) 

(11) 
(12) 



The following theorem can be proved by a direct calcu- 
lation (for a simple proof see 



Theorem 1. 



id a (m\ = m\(v[i]-xj), 

-id s m = m(v[i]-\j). 



(13) 
(14) 



Remarks: (a) The above transformation is called "bi- 
nary" because it can be shown to be a composition of two 
mutually conjugated "elementary" Darboux-type trans- 
formations introduced and analyzed in In this re- 
spect it is similar to the standard binary Darboux trans- 
formation discussed in 

(b) The form of the transformation resembles the ex- 
pressions occuring in the so-called dressing method intro- 
duced by Zakharov in the context of the Riemann-Hilbcrt 
factorization problem |j] |llj . The difference between the 
binary transformation and the dressing method lies es- 
sentially in the form of P and restrictions imposed on 
the parameters p and v. In the original Zakharov con- 
struction the parameters were related to poles of complex 
functions and the whole approach relied on analyticity 
properties of objects under consideration. This complex- 
analytic framework made it difficult to work with infinite- 
dimensional systems. The binary Darboux technique is 
purely algebraic and therefore the above restrictions no 
longer apply. The parameters are arbitrary and an ex- 
tension to infinite-dimensional systems is natural, cf. the 
example analyzed in (5) where J = H was a Hamiltonian 
of a quantum-mechanical harmonic oscillator. 



III. NONLINEAR DENSITY MATRIX 
EQUATIONS 

In this section we discuss a class of nonlinear equations 
that can be solved by the binary Darboux technique de- 
scribed in the previous section. Their common feature is 



the existence of a Darboux-covariant Lax pair. The class 
is parametrized by a natural number n, and for n = 1 
reduces to EAvNE analyzed in ||. The EAvNE for pro- 
jectors reduces to the pure-state vNE which, on the other 
hand, is equivalent to the Schrodinger equation (SE). For 
n > 1 the equations we discus in this section are equiva- 
lent, for pure states, to a class of nonlinear Schrodinger 
equations (NSE). 

Let A be a time-independent self-adjoint operator. 
Consider the following class of nonlinear vNE's 



ip = [H(p),p] 



J2[A n - k P A k ,p] 

k=0 



k=0 



(15) 
(16) 



where the dot denotes the time derivative. For n = 1 and 
A = H one finds 



H(p) = Hp + pH 
and the nonlinear vNE becomes the EAvNE 
ip=[Hp + pH,p] = [H,p 2 ]. 



(17) 



(18) 



Let us stress that the problem we consider is generically 
infinite-dimensional and p is positive, bounded, trace- 
class and Hermitian. These constraints are Darboux- 
covariant, as will become clear later, if the binary v = pZ 
transformation is considered. 

The Lax representation of (|l^) is given by the pair of 
linear equations 



z M = (p- v a )\<p)i 

n 

i\ip) = (Y / A n - k P A k -pA n+1 )\v), 



(19) 
(20) 



k=0 



where z^, p, € C. In order to use the binary technique 
we have to consider the additional two conjugated Lax 
pairs with parameters A, Z\ , v, z v : 



= {^\{p-\A), 

n 

= {il)\(X2A n - k pA k - \A n+1 } 

k=0 

z v{x\ = (x\{p - vA), 

n 

-i(x\ = (x\(J2 An ~ k P Ak - vAn+1 



(21) 
(22) 
(23) 
(24) 



k=0 



In what follows the p- and ;/-pairs (|19[)— ( 20 ) and (|2ji 
( p4| ) will be used to define the binary transformation of 
the conjugated A-pair (|2l|)-(p2]). 

Lemma 1. The Lax pair (pl|) (p2|) is covariant under 
the binary Darboux transformation. 

Proof: The Lax pair can be regarded as a particular case 
of the zero-curvature representation 
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-i(ip\ 



(ip\(p-XA) =: (^K^-AJi), (25) 

n 



fe=0 



=: (V'K^-AJs), 



(26) 



where the prime denotes a derivative with respect to some 
auxiliary parameter r, and // = 0. Eqs. (P5|), ( p6|) are 
examples of general Zakharov-Shabat problems. The bi- 
nary Darboux transformation implies 



where 
Vi[l] = 

V*[l] = 



i{m'\ = (m\{vi[i]-xji) 
-i(m\ = (^[i]|(f 2 [i]-aj 2 ) 



p + (p-u)[P,A] =:p[l] 
V 2 + {p-v)[P,J 2 ] 



(27) 
(28) 



(29) 
(30) 
(31) 



in+ll 



Y A n - k P A k + (p - v)[P, A n 

k=0 

n n 

Y A n ' k pA k + {p-v)J2 A n - k [P, A]A k 

k=0 k=0 
n 

J2A n - k p[l]A k . (32) 



fe=0 



The r-dependence of the solutions of the zero-curvature 
equations is given by 



|^,T))=e-^>(t,T = 0)) 
( x (t, T )\=e i ^(x(t,r = 0)\ 

m,T)\ = e ^m,T = o)\ 

so that P is r-independent and, hence, 

^[l](t,r)|=e te ^(^[l](t,r = 0)|. 
Therefore 

-i(m'\ = (m\m-*A)= Zx (m 

71 

- t (m\ = (m\(Y, An ~ k pW Ak - XA 

k=0 



n+1 



(33) 
(34) 
(35) 



(36) 

(37) 
(38) 



which was to be proved. □ 

Remarks: (a) The above Lemma implies that Eqs. Jl6| ) 
are invariant under the action of the binary Darboux 
transformation, since the stationarity conditions imply 
P' = and p[l}' = if p' = 0. 

(b) For v — p p = p^ implies p[l] = 

(c) Trp= Trp[l]. 

(d) For pure states, i.e. p = |*)(*|, (*|*) = 1, 
Eqs. ( [l6| ) are equivalent to a class of nonlinear SE. In- 
deed, in this case 

pA k p= |*)(*|A fc |*)(*| = Tr(pA k )p 



and 



ip = P A k p] = Y, Tr ( P A k )[A n - k , p}. (39) 

k=0 k=0 

The solutions of (|3|) are 

since Ti (pA k ) are time- independent. The same dynam- 
ics could be obtained from the NSE 



i*) = Y^ Ak \^ An ~ k \^- 



(40) 



k=0 



Let us recall that the case n = 1, A = H is equivalent 
to the EAvNE which for pure states gives the ordinary 
vNE. This property is seen also in ( [40| ) since then the 
corresponding SE is the ordinary linear one. 



IV. BINARY DARBOUX TRANSFORMATION 
AS A SIMILARITY TRANSFORMATION 

Spectra of all the solutions of the EAvNE discussed 
in Q were invariant under the binary transformation. 
Below we show that the result is not a pecularity of the 
EAvNE but holds for all operator equations that are com- 
patibility conditions for a Lax pair of ZS equations. In 
a more general situation where instead of a Lax pair we 
have a zero-curvature pair, the result does not have to 
hold. For solutions that are not self-adjoint the result 
extends separately to both the left and the right spectra. 

Consider the following three general zero-curvature 
pairs 



i-W) = {p- P A )W) = 

= (V(p)-pJ)\if), 

-i{x'\ = (x\(p ~ vA) = z v {x\ 
-i{x\ = (x\{V(p)-uJ), 
-i(1/\ = (1>\(j>-\A) 
-m = (iP\(V(p)-\J), 



(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



The only assumption we make about p, J, A, and V(p) 
is the covariance of (|H|)-(^) under the binary Darboux 
transformation constructed with the help of (x\ and \ip). 
Let us note that the stationarity with respect to r is 
assumed only for the first two pairs. 

Theorem 2. Under the above assumptions the binary 
Darboux transformation p p[l] is a similarity trans- 
formation, p[l] = TpT^ 1 , where 



(47) 



Proof: By definition p[l] = p + (p - v)[P,A}. Eqs. (|4lj), 
@ imply 
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Z V P 

from which it follows that 



(P - »A)P 
P{p - vA) 



P{ P - pA)P = { P - pA)P 
P(p-vA)P = P(p-vA). 



(48) 
(49) 



(50) 
(51) 



Px(t) 



-i(n+l)XA n t 



{p(t)+X)e 



i(n+l)XA"t 



(54) 



Multiplying ( |5C| ) by v, (|5lJ) by p, and subtracting the 
resulting equations we get 



V — LI 1 1 

[P, A] = -PpP - -pP + -Pp. 

pv p V 



(52) 



Inserting this expression into the definition of p[l] we 
obtain 



P[l] 



l + ^p)p(l + ^P 



(53) 



To see that the operators occuring at both sides of p are 
inverses of each other it is sufficient to use P 2 — P. For 
any z £ C and any P satisfying P 2 = P one finds 

exp(zP) = 1 - P + e z P 



and, in particular, 
□ 



1 - P + -P = T. 

v 



Remarks: (a) For Hcrmitian P and v = p the operator 
T is unitary. If in addition p is imaginary then T 
I 21'. 

(Jo) The assumptions we have made can be weakened 
since the proof is based on the equalities (50) and j5l| ) 
which are weaker than (^d]), (|43|). 

(c) The proof extends without any modification to the 
operator version of the binary Darboux transformation, 
where we take the operators (p and x instead of the 
Hilbert-space vectors \ip) and (x\, and P — (^(px^p)^ 1 ^- 

(d) The result applies to all the equations described in 
the previous section. As a consequence all their solutions 
p[l] obtained from a density matrix p with v — p are also 
density matrices. 



V. SPECTRUM SHIFTING AND RESCALING 

It is often convenient to start with a solution that is 
not a density matrix. For example, one can start with a 
p which is not positive, or whose trace is not equal to 1. 
To turn it into a density matrix one first has to add an 
appropriate positive operator to p and then renormalize 
the solution to get Trp = 1. Below we give examples of 
such transformations. 

Take an operator X satisfying [A, A] — [X,p] = 0. 
Such an X always exists (e.g. X = Al, A e R). Assume 
p(t) is a solution of (|l6|). Then 



also satisfies (|l6|). In particular, for X — Al the spec- 
trum of px is shifted by A with respect to this of p. 
For n = 1 one obtains the gauge transformation of the 
EAvNE discussed in || . This trick is very useful in prac- 
tical computations of finite dimensional cases and typi- 
cally leads to positive but non-normalized solutions. 

To change the normalization one uses the rescaling 
symmetry of (|l6|): 



p(t)~ PY (t) = Yp(Yt), Y e R. 



(55) 



As an application of the above two transformations con- 
sider the problem of stationary seed solutions p. The so- 
lutions discussed in were starting with non-stationary 
p's, since those satisfying [p, A] = lead to the projector 
P commuting with both A and p, and the binary trans- 
formation is trivial. Still, there exists another class of 
stationary solutions of (|16|), obtained if Ap = —pA. Now 
the projector P will not, in general, commute with p and 
A, and the binary transformation may be nontrivial. 

Assume [P,p] ^ 0. For n = 2m one has X)fc=o( — V) 
1 and 



#) = (A>-M" +1 )|V> = V4>) 
whose formal solution is 

\<p(t))=e- i *» Ant M0)). 
For odd n we get ^2^, =0 (— 1) = and 
i\<p) = -pA n+1 \ v ) 



implying 



e^" +lt |^(0)). 



(56) 



(57) 



(58) 



(59) 



The fact that p and z M are complex results, in general, in 
a nontrivial time evolution of P arising from the nontriv- 
ial contribution from the time dependence of (ip(t)\tp(t)) . 
Let us note here that the "phase transitions" found for 
the Euler-Arnold-von Neumann equation had a similar 
origin although the "seed solutions" p satisfied a differ- 
ent condition. An exceptional situation occurs if the so- 
lution of the Lax pair turns out to be an eigenstate of 
A 2m (for n = 2m) or A 2( - m+1 ^ (for n = 2m + 1). Let 
us also note that p that anticommutes with A does com- 
mute with A 2 and therefore also with the generators of 
the time evolution given by (|57|) and (|59|). This is an 
important property allowing to look for an eigenstate of 
p — pA at t = 0. 

For p anticommuting with A one finds Tr pA = 0. This 
implies that p and A cannot be simultaneously positive. 
In a finite dimensional case one can easily shift the spec- 
trum of p by a number and then renormalize the resulting 
solution by means of the above two transformations, and 
in this way use p as an intermediate step towards a den- 
sity matrix p[l] with positive A. An analysis of explicit 
solutions obtained by this technique will be presented in 
a forthcoming paper [O . 
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VI. FURTHER PERSPECTIVES 

The central result of this Letter is the formula ( ^3| ) . For 
\i 7^ v the operator T is nonunitary. A possibility that has 
not been explored as yet is to look for Hermitian p[l] but 
starting with a non-Hermitian p and nonunitary T. In 
still more general perspective, it is clear that the method 
we have developed works for a class of nonlinear operator 
equations of the Heisenberg form iF = [H(F),F] and, 
hence, may find applications in quantum optics and field 
theory. The technique is well suited for solving infinite- 
dimensional problems, but a practical complication is 
associated with constraints typical of density matrices. 
Hermiticity is an easy issue here, but real difficulties oc- 
cur if one looks for trace-class and positive solutions. 

Another open question is how to solve vNE's of the 
general form 

ip=[H,f(p)}. (60) 

There are reasons to investigate this type of nonlincaritics 
since there exist links of f(p) = p q , q > 0, to nonextensive 
g-statistics (cf. jl3| and references therein). At the mo- 
ment no binary-Darboux covariant Lax pairs are known 
for g^2. Actually, looking more closely at the struc- 
ture of Zakharov-Shabat equations and properties of the 
binary transformation one can conjecture that this par- 
ticular Darboux transformation can work only for vNE's 
with Hamiltonians linear in p. For example, we are aware 
of a Lax pair leading to 

lp = [H, p 3 } = [p 2 H + pHp + Hp 2 , p] = [H{p),p] (61) 

but it is not covariant under the binary transformation: 
The reason is the quadratic dependence on p of the 
Hamiltonian H(p). On the other hand, all /(/?) nonlin- 
earities with polynomial / have similar properties when 
looked at from a Lie-Poisson perspective (say, the same 
Casimirs and other invariants) and there is basically no 
reason why these equations should not be solvable by a 
Darboux-type technique. It is possible that such trans- 
formations should be looked for between those obtained 
by a composition of three or more elementary transfor- 
mations discussed in @||. 

A new class of constraints, which have not been consid- 
ered so far, are those related to commutation relations. 
To give an example, in quantum optics the task is of- 
ten reduced to solving nonlinear Heisenberg equations for 
creation and annihilation operators. Assuming the prob- 
lem is treated by a Darboux technique, the constraint 
will be of the form [a[l],a[l]t] = 1 and the fact that we 
are dealing with a similarity transformation (^3|) may be 
of crucial importance. 

All these problems are a subject of current study. 
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